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IX, Abstract 

In this paper, we study the number of different interference alignment (IA) solutions that exists for 
the if -user multiple-input multiple-output (MIMO) interference channels with constant coefficients, when 
c/3 the alignment is performed via beamforming and without symbol extensions. When counting the number 

of IA solutions for a given scenario, the most interesting case arises when the number of equations in 
the polynomial system matches the number of variables and the system is feasible. In this situation, 

> 

the number of IA solutions is finite and constant for any channel realization out of a zero-measure set 

ON 



and, as we prove in the paper, is given by an integral formula that can be numerically approximated 
using Monte Carlo integration methods. More precisely, the number of alignment solutions is the scaled 



average over a subset of the solution variety (formed by all triplets of channels, precoders and decoders 
satisfying the IA polynomial equations) of the determinant of a certain Hermitian matrix related to the 
geometry of the problem. Interestingly, while the value of this determinant at an arbitrary point can be 
^ used to check the feasibility of the IA problem, the average of the determinant (properly scaled) gives 

us the number of solutions. Our results can be applied to arbitrary interference MIMO networks, with 
any number of users, antennas and streams per user. 
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I. Introduction 

Interference alignment (IA) has received a lot of attention in recent years as a key technique to achieve 
the maximum degrees of freedom (DoF) of wireless networks in the presence of interference. Originally 
proposed in (T], (2J, the basic idea of IA consists of designing the transmitted signals in such a way that 
the interference at each receiver falls within a lower-dimensional subspace, therefore leaving a subspace 
free of interference for the desired signal pj. This idea has been applied in different forms (e.g., ergodic 
interference alignment [4], signal space alignment |TJ, or signal scale alignment J5J, (6J), and adapted 
to various wireless networks such as interference networks |1], X channels (2j, downlink broadcast 
channels in cellular communications [7] and, more recently, to two-hop relay-aided networks in the form 
of interference neutralization (8j. 

In this paper we consider the linear IA problem (i.e., signal space alignment by means of linear 
beamforming) for the i^-user multiple-input multiple-output (MIMO) interference channel with constant 
channel coefficients. Moreover, the MIMO channels are considered to be generic, without any particular 
structure, which happens for instance when the channel matrices have independent entries drawn from a 
continuous distribution. This setup has been the preferred option for recent experimental studies on IA 

®, (To), (nj. 

The feasibility of linear IA for MIMO interference networks, which amounts to study the solvability 
of a set of polynomial equations, has been an active research topic during the last years p2j , (T3J, 
(14) , p3| , [ [To*} . Combining algebraic geometry tools with differential topology ones, it has been recently 
proved in fTf] , that an IA problem with any number of users, antennas and streams per user, is feasible 
iff the linear mapping given by the projection from the tangent space of V (the solution variety, whose 
elements are the triplets formed by the channels, decoders and precoders satisfying the IA equations) to 
the tangent space of H. (the complex space of MIMO interference channels) at some element of V is 
surjective. Note that this implies in particular that the dimension of V must be larger than or equal to 
the dimension of T~L. Exploiting this result, a general IA feasibility test with polynomial complexity has 
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also been proposed in |17|, [18|. This test reduces to check whether the determinant of a given square 
Hermitian matrix is zero (meaning infeasible almost surely) or not (feasible). 

In this paper we study the problem of how many different alignment solutions exist for a feasible IA 
problem. While the number of solutions is known for some particular cases, a general result is not available 
yet. For instance, it can be trivially shown that feasible systems for which the algebraic dimension of the 
solution variety is larger than that of the input space have an infinite number of alignment solutions. In 
plain words, these are MIMO interference networks for which the number of variables is larger than the 
number of equations of the polynomial system. These scenarios typically represent cases where not all 
available DoF are achieved (for instance, we might have more antennas than strictly needed to achieve 
a certain DoF tuple) and therefore they do not receive further consideration in this paper. Much more 
interesting and challenging is the case where the dimensions of V and H. are exactly the same (identical 
number of variables and equations) and the problem is feasible, because in this situation the number 
of IA solutions is finite and constant out of a zero measure set of H. as proved in fTTJ . Following the 



nomenclature recently introduced in |19| we refer to these systems as tightly feasible, stressing the fact 
that removing a single antenna from the network turns the IA problem infeasible. 

For tightly feasible single-beam (i.e., when all users wish to transmit d = 1 stream of data) MIMO 
networks, and elaborating on classic results from algebraic geometry, it was shown in [12] that the 
number of alignment solutions coincides with the mixed volume of the Newton polytopes that support 
each equation of the polynomial system. Although this solves theoretically the problem for single-beam 
networks, in practice the computation of the mixed volume of a set of IA equations using the available 



software tools |20| can be very demanding, therefore only a few cases have been solved so far. For 
single-beam networks, some upper bounds on the number of solutions using Bezout's Theorem have also 



been proposed in [ 12], [21 ]. For multi-beam scenarios, however, the genericity of the polynomials system 
of equations is lost and it is not possible to resort to mixed volume calculations to find the number of 
solutions. Furthermore, the existing bounds in multi-beam cases are very loose. 

The main contribution of this paper is an integral formula for the number of IA solutions for arbitrary, 
tightly feasible, networks. More specifically, we prove that while the feasibility problem is solved by 
checking the determinant of a certain Hermitian matrix, the number of IA solutions is given by the 
integral of the same determinant over a subset of the solution variety scaled by an appropriate constant. 
Although the integral, in general, is hard to compute analytically, it can be easily estimated using Monte 
Carlo integration. To speed up the convergence of the Monte Carlo integration method, we specialize 
the general integral formula for square symmetric multi-beam cases (i.e., equal number of transmit and 
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receive antennas and equal number of streams per user). 

In addition to having a theoretical interest, the results proved in this work might also have some 
practical implications. For instance, to find scaling laws for the number of solutions with respect to the 
number of users could have interest to analyze the asymptotic performance of linear IA, as discussed 



in [21 1. Also, for moderate-size networks for which the total number of solutions is not very high, the 
results of this paper also open the possibility to provide a systematic way to compute all (or practically 
all) interference alignment solutions for a channel realization. This idea is also briefly explored in the 
paper. 

The rest of the paper is organized as follows. In Section [III the system model and the IA feasibility 



problem are briefly reviewed, paying special attention to the feasibility test in |17|. The main results of 
the paper are presented in Section [Illj where two integral formulas for the number of IA solutions are 
provided, one valid for arbitrary networks and the other for symmetric multi-beam scenarios (all users 
with the same number of antennas at both sides of the link and transmitting the same number of streams). 
Although these integrals, in general, cannot be computed in closed form, they can easily be estimated 
using Monte Carlo integration. A short review on Riemmanian manifolds and other mathematical results 



that will also be used during the derivations as well as the proofs of the main theorems in Section III 
are relegated to appendices. 

II. System model and background material 

In this section we describe the system model considered in the paper, introduce the notation, define 
the main algebraic sets used throughout the paper, and briefly review the feasibility conditions of linear 
IA problems for arbitrary wireless networks. 



A. Linear IA 

We consider the K-user MIMO interference channel with transmitter k having > 1 antennas and 
receiver k having > 1 antennas. Each user k wishes to send df. > streams or messages. We 
adhere to the notation used in p2| and denote this (fully connected) asymmetric interference channel 
as Y\k=i x Nk, dk) = (Mi x N\,d\) ■ ■ ■ (Mk X Nk, dx)- The symmetric case in which all users 
transmit d streams and are equipped with M transmit and N receive antennas is denoted as (M x N,d) K . 
In the square symmetric case all users have the same number of antennas M = N. 

The MIMO channel from transmitter / to receiver k is denoted as H^i and assumed to be flat-fading 
and constant over time. Each fljy is an x M\ complex matrix with independent entries drawn from 
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a continuous distribution. We denote the set of users as T = { 1 , . . . , K} and the set of interfering links 
as $> = {(k, I) G T x T : k 7^ I}. Also, denotes the cardinality of <I>, that is the number of elements 
in the finite set <I>. In this paper we focus on fully connected interference channels and, consequently, 

User j encodes its message using an Mj x dj precoding matrix Vj and the received signal is given by 

Uj = HjjVjXj + ^ H ji V i x i + n j, l<j<K (1) 

where xj is the dj x 1 transmitted signal and nj is the zero mean unit variance circularly symmetric 
additive white Gaussian noise vector. The first term in ([j]) is the desired signal, while the second term 
represents the interference space. The receiver j applies a linear decoder Uj of dimensions Nj x dj, i.e., 

Ujyj = UjHjjVjX, + £ UfHjiViXi + Ujrij , 1 < j < K, (2) 

where superscript T denotes transpose. 

The interference alignment (IA) problem is to find the decoders and precoders, Vj and Uj, in such a 
way that the interfering signals at each receiver fall into a reduced-dimensional subspace and the receivers 
can then extract the projection of the desired signal that lies in the interference-free subspace. To this 
end it is required that the polynomial equations 

U^HmV^O, (M)€$, (3) 

are satisfied, while the signal subspace for each user must be linearly independent of the interference 
subspace and must have dimension d^, that is 

i^Y(UlH kk V k ) = d k , V k € T. (4) 

B. Feasibility of IA: a brief review 

The IA feasibility problem amounts to study the relationship between dj , Mj ,Nj,K such that the linear 
alignment problem is feasible. If the problem is feasible, the tuple {d\ , . . . , dx) defines the degrees 
of freedom (DoF) of the system, that is the maximum number of independent data streams that can 
be transmitted without interference in the channel. The IA feasibility problem and the closely related 
problem of finding the maximum DoF of a given network have attracted a lot of research over the last 
years. For instance, the DoF for the 2-user and, under some conditions, for the symmetric K-user MIMO 



interference channel have been found in [22] and |23J, respectively. In this work we make the following 
assumptions: 

1 < d k < N k , V k e T, 1 < di < M h V I e T, (5) 
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and 

d k + di<N k + M h V(M)€$, (6) 

which are necessary conditions for feasibility derived, respectively, for point-to-point MIMO links and 
for the 2-user MIMO channel. 



The IA feasibility problem has also been deeply investigated in [12|-[16|. In the following we make a 
short review of the main feasibility result presented in (T7J, fT8[ , which forms the starting point of this 
work. 

We start by describing the three main algebraic sets involved in the feasibility problem. 
• Input space formed by the MIMO matrices, which is formally defined as 

H= [J M NkXMl (C) (7) 
(fc,l)e* 

where Y\ holds for Cartesian product, and M-N h xMi( l C) is the set of N k X Mi complex matrices. 



Note that in [ 17 1, [ 18 1, we let % be the product of projective spaces instead of the product of affine 
spaces. The use of affine spaces is more convenient for the purposes of root counting. 
Output space of precoders and decoders (i.e., the set where the possible outputs exist) 



Vfeex / vex / 



(8) 



where G a & is the Grassmannian formed by the linear subspaces of (complex) dimension a in C b . 
The solution variety, which is given by 



V = {(H, U,V)eHxS: © holds} (9) 

where H is the collection of all matrices H k i and, similarly, U and V denote the set of U k and Vi, 
respectively. The set V is given by certain polynomial equations, linear in each of the H k i,U k ,Vi 
and therefore is an algebraic subvariety of the product space % x S. Let us remind here that 
the IA equations given by ([3]) hold or do not hold independently of the particular chosen affine 
representatives of U, V. 

Once the main algebraic sets have been defined, it is interesting to consider the following diagram 

V 

fl / \ 7T 2 (10) 

u s 
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where the sets and the main projections involved in the feasibility problem are depicted. Note that, given 
H G U, the set vrf 1 (H) is a copy of the set of U, V such that ^ holds, that is the solution set of the 
linear interference alignment problem. On the other hand, given (17, V) G S, the set tt^{U,V) is a copy 
of the set of H G Ti such that ((3) holds. 

The feasibility question can then be restated as, is TT^iH) / %for a generic H? The question was 
solved in p7| , basically stating that the problem is feasible if and only if two conditions are fulfilled: 

1) The algebraic dimension of V must be larger than or equal to the dimension of T~L, i.e., 



Y,d k (N k + M k -2d k ) 1 - Yl d kdi\>0- 

aeT / \(k,i)e-s> J 



(11) 



In other words this condition means that, for the problem of polynomial equations to have a solution, 
the number of variables must be larger than or equal to the number of equations. This condition 



was already established in |12|, hereby classifying interference channels as proper (s > 0) or 
improper (s < 0). More recently, in p3| it was rigorously proved that improper systems are always 
infeasible. 



2) For some element (H, U, V) G V, the linear mapping 



({Uk}ker,{Vi}ier) -> {^H^Vi + U^HkM} 



(12) 



(fe,0e* 

is surjective, i.e., it has maximal rank equal to Yl(k 0e* dkdi- This condition amounts to saying that 
the projection from the tangent plane at an arbitrary point of the solution variety to the tangent plane 
of the input space must be surjective: that is, one tangent plane must cover the other. Moreover, in 



this case, the mapping ( 12 1 is surjective for almost every (H, U, V) G V. 



III. The number of solutions of feasible IA problems 

A. Preliminaries 



As it was shown in |17|, [18], the surjectivity of the mapping 9 in (12i can easily checked by a 
polynomial-complexity test that can be applied to arbitrary TT-user MIMO interference networks. The test 
basically consists of two main steps: i) to find and arbitrary point in the solution variety and ii) to check 
the rank of a matrix constructed from that point. To find an arbitrary point in the solution variety, in |T7j 
we generated a set of random precoders and decoders, and then solved a linear underdetermined problem 
to get a set of channel matrices satisfying the IA equations ([3]): this was called inverse I A problem in fT7| , 
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|18|. In this paper we choose an even simpler (trivial) solution satisfying the I A equations. Specifically, 
we take structured matrices given by 

Od k xd t A k i 
Bki Cki 

with precoders and decoders given by 

Id, 



H, 



hi 



(13) 



Vi 



(Mi-di)xdi 



(N k -d k )xd k 



(14) 



which trivially satisfy U^H^Vi = and therefore belong to the solution variety. We claim that essentially 
all the useful information about V can be obtained from the subset of V consisting on triples (Hki, Uk,Vi) 



of the form ( 13 1 and ( 14 1. The reason is that given any other element (Hki, Uk, V{) G V, one can easily 



find sets of orthogonal matrices Pk and Q\ satisfying 

U k = PkUk, Vi = QiVi, 



and 



UlHuVi = U'£ (Pi) 1 H k iQ*Vi = 0, 



tT 



T fV 



where the superscript * denotes Hermitian. That is, the transformed channels Hki = (Pf,) T HkiQ*i nave 



the form ( [13] ), and the transformed precoders V\ and decoders Uk have the form ( [T4| ). Thus, we have just 
described an isometry which sends (Hki, Uk, Vi) to (Hki, Uk, Vi). The situation is thus similar to that of 
a torus: every point can be sent to some predefined vertical circle through a rotation, thus the torus is 
essentially understood by "moving" a circumference and keeping track of the visited places. The same 



way, V can be thought of as moving the set of triples of the form (13) and (14), and keeping track of 
the visited places. Technically, V is the orbit of the set of triples of the form (13) and ( [T4] ) under the 
isometric action of a product of unitary groups. 

In summary, the main idea is that, for the purpose of checking feasibility or counting solutions, we 
can replace the set of arbitrary complex matrices % by the set of structured matrices 

Bki Cki 




7T. 



-1 



h k 



\Q(N k -d k )xd k J fceT K^iMi-d^xdJ leT/ 



Id, 



The mapping 6 in ( 12 1 has a more simple form for triples of the form ( 13 1 and ( 14 1, and can be replaced 



by a new mapping ^ defined as 

* : {Uker M (N k -d k) xd k (C)) x (UierM {Ml -dOxdX c )) 
({Ukher, {Viher) 



U( k ,i)^ M d k xdAC) 

k 



VTBu + A M Vi 



(15) 



k.i 
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We will be interested in the function det^^*), which depends on the channel realization H only through 
the blocks A k \ and B k \. The vectorization of the mapping ([13]) reveals that \t is composed of two main 



kinds of blocks, and ^7 , i.e. 

T ; 



(A) 



vec(U£ B kl + A M Vi) = (A kl ® I dk )K (Nk _ dk)jdk ^ (U k ) + (J d , Bfr) vec(Vj), 
" (A) 



(16) 



' kl ~kl 
where ® denotes Kronecker product and K m ^ n is the ran x mn commutation matrix which is defined 

as the matrix that transforms the vectorized form of an m x n matrix into the vectorized form of its 

transpose. Block Vl/j^ has dimensions did k x di{M\ — d{), whereas block is d\d k x d k (N k — d k ). 

For a given tuple (k, I), ^ k f^ and ^If are placed in the row partition that corresponds to the interfering 

link indicated by the tuple (k, I). V^jy is placed in the l + K-th column partition, whereas occupies 

the k-th column partition. The rest of blocks are occupied by null matrices. The dimensions of ^ are 

therefore J2k^i ^kd-i x Y^j=iiMj + ~ %dj)dj. In the particular case of s = 0, \& is a square matrix 

of size ^2 k -£id k di. 

Notice that ^ has the same structure as the incidence matrix of the network connectivity graph. Taking 
the 3-user interference channel as an example, ^ is constructed as follows 



(A) 
12 

(A) 
13 












,{A) 
21 

,{A) 
■23 
















(B) 
21 





(A) y(B) 



31 

(A) 
32 



31 





(B) 

12 






(B) 
32 



(B) 
13 



(B) 
23 






where the blocks Vpjy and are given by ( fT6] ). 



(-4) 



B. Mam results 

We use the following notation: given a Riemannian manifold X with total finite volume denoted as 
Vol(X) (the volume of the manifolds used in this paper are reviewed in Appendix |A|), let 

/ f(x) dx = 1 / f(x) dx 

be the average value of a integrable (or measurable and nonnegative) function / : X— >M.. Fix dj,Mj, Nj 
and <3? satisfying ([5]) and ([6]) and let s > be defined as in (11). The main results of the paper are 
Theorems 1 , 2 and 3 below, which give integral expressions for the number of IA solutions when s = 
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and the system is tightly feasible: this number is denoted as ^(tt^ 1 (Hq)), which is the same for all 
channel realizations out of some zero-measure set. 

Theorem 1: Assume that s = 0, and let Tie C % be any open set such that the following holds: if 
77 = (Hki) £ H e and P^, Qk, 1 < k < K are unitary matrices of respective sizes TVjt, Mfc, then 

(We may just say that % e is invariant under unitary transformations). Then, for every Hq £ % out of 
some zero-measure set, we have: 

^ 1 (H )) = C [ det(W) d£7, (17) 

where 

y z(?7 e )' 

with 5 being the output space (Cartesian product of Grasmannians) in Eq. ([8]). 

Proof: See Appendix [B] ■ 
If we take H e to be the set 

{(H M ) : \\H M \\ F E (l- ej l + e )} 

and we let e i— > we get: 

Theorem 2: For a tightly feasible (s = 0) fully connected interference channel, and for every Hq £ % 
out of some zero-measure set, we have: 

^ 1 (H )) = C-f det(***) dH, 

JHeHi,\\H M \\ F =l 

where 

r(jv t M,) 



c- n 

(fc,0e* 



r(iv fc M, - 
r(2)---r(4)-r(2)---r(iV fc -4) 



2 V l'(2)...L-(A-,) 



n 

Proof: See Appendix [C] 



r(2)---r(d;)-r(2)---r(M;-di) 
r(2)---r(Af { ) 



Remark 1: As proved in |T7J (see also |T8J), if the system is infeasible then det^^*) = for every 
choice of 77, U, V and hence Theorem [T] still holds. Moreover, if the system is feasible and s > then 
there is a continuous of solutions for almost every 77^ and hence it is meaningless to count them (the 
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value of the integrals in our theorems is not related to the number of solutions in that case). Note also 
that the equality of Theorem [1] holds for every unitarily invariant open set T-L e , which in particular implies 
that the right-hand side of ( [17) ) has the same value for all such 7-L € (recall that we proved in [17] that 
almost all channel realizations in % have the same number of solutions). 

Theorem [2] can be used to approximate the number of solutions of a given MIMO system using 



Montecarlo integration (see Section III-C below). However, the convergence of the integral is quite slow 
in general. In the square symmetric case when all the dp. and all the and are equal \/k and greater 
than 2d, that is N = M > 2d, which holds automatically when s = and K > 3; we can write another 
integral which has faster convergence in practice: 

Theorem 3: Let us consider a tightly feasible (s = 0) square interference channel (N^ = = N 
and dk = d, \/k). Assuming additionally that K > 3, then for every Hq £ % out of some zero-measure 
set, we have: 

= i v 2 2lvwf v Y*W)/ det(W) dH, 

\Vol(U N )Vol(U N -2d) J hA> M ,B hl )ati N _ md 

where ^ is again defined by ( [13] ) and the input space of MIMO channels where we have to integrate are 

now 



ki 



B kl Q(N-d)x(N-d). 



whose blocks, A^i and B^, are matrices in the complex Stiefel manifold, denoted as U(N-d)xd> an d 
formed by all orthonormal (f-dimensional vectors in C^ N ~ d \ On the other hand, U a denotes the unitary 
group of dimension a, whose volumen can be found in Appendix [Aj 

Remark 2: If the problem is fully connected, the value of the constant preceding the integral in Theorem 

[3] is (using that 2N - dK - d = if s = 0): 

J2 \ K(K-l) 

Vol{U N )Vol{U N . 2d ) ) 1 ' 

■ ■ t(n) \ K ( K ~ l ) ( r(2)---r(d) 




T(N - 2d + 1) • • • T(N -d)J \T(N - d + 1) • • • T(N), 

Additionally, if N = 2d (which implies K = 3) then this constant is exactly equal to 1. 

Proof: See Appendix [D] ■ 
In the next section we discuss how the results in Theorems 1 and 2 can be used to get approximations 
to the number of IA solutions for a given interference network. 
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C. Estimating the number of solutions by Monte Carlo integration 

The integrals in Theorems [2] and [3] are too difficult to be computed analytically, but one can certainly 
try to compute them approximately using Monte Carlo integration. Our main reference here is [24, Sec. 
5]. The Crude Monte Carlo method for computing the average 

f f(x)dx 

JxeX 

of a function / defined on a finite- volume manifold X consists just in choosing many points at random, 
say xi, . . . ,x n for n » 1, uniformly distributed in X, and approximating 



r 1 n 

/ f(x)dx^E n = -Y j f(, 
Jxex n 



(18) 



The most reasonable way to implement this in a computer program is to write down an iteration that 
computes E\, E2, -E3, . . .. The unique point to be decided is how many such Xj we must choose to get a 
reasonable approximation of the integral. A usual tool for measuring that is the standard deviation, that 
can be approximated by 



I n — 1 

\ j=i 



S n = I > .(/(■'•,) - t„)~ ) • (19) 

If we stop the iteration when 




E n 



then, with a probability of 0.95 on the set of random sequences of n terms, the relative error satisfies 

\I xeX f( X ) dx - E n\ 



< 2e. 



\E n \ 

For example, if we stop the iteration when 

^ < 0.05 

then we can expect to be making an error of about 10 percent in our calculation of f xl - x f(x) dx. 
The whole procedure for a general system is illustrated in Algorithm [T] which follows Theorem [2] Its 
particularization to square systems is shown in Algorithm [2] and follows Theorem [3] 

D. The single-beam case 

Although the results of Theorems 1, 2 and 3 are general and can be applied to arbitrary systems, for the 
particular case of single-beam MIMO networks (dk = 1, k € T) it is possible to develop specific, much 
more efficient, techniques to count the exact number of alignment solutions. This subsection is devoted 
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Algorithm 1: Computing the number of IA solutions for general scenarios Y\k=i iMk x ^k-,dk) 
Input: Relative error, e; number of antennas, {Mi} and {iVfc}; streams, d; and users, K. 

Output: Approximate number of I A solutions, E n . 

n = 1 

repeat 

Generate a set of i.i.d. matrices {A^i}, {B^} and {C^}. 
Build channel matrices {H^i} according to ( p"3] ). 
Normalize every channel matrix Hki such that = 1. 

Build matrix ^ according to ( [To] ). 

Compute D n = C det(^^f*) where C is taken from Theorem j^J 

Calculate E n and S n according to (_18i and ([19]), respectively, where f(xj) is now Dj. 



n = n + 1. 
until |f < e 



Algorithm 2: Computing the number of IA solutions for symmetric square scenarios (N x N, d) 
Input: Relative error, e; number of antennas, N; streams, d; and users, K. 

Output: Approximate number of I A solutions, E n . 

n = 1 

repeat 

Generate a set of i.i.d. (N — d) x d matrices {^4.t;} and {B^} in the Stiefel manifold. 
Build matrix ^ according to ( fT6[ ). 

Compute D n = Cdet(*I'^'*) where C is taken from Theorem j^J 



K 



Calculate E n and S n according to ( fT8| ) and ([19]), respectively, where f(xj) is now Dj. 

n = n + 1. 
until |f- < e 



to this particular case. First, we should mention that, from a theoretical point of view, the single-beam 



case was solved in [12|, where it was shown that the number of IA solutions for single-beam feasible 
systems coincides with the mixed volume of the Newton polytopes that support each equation of the 
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system^ However, from a practical point of view, the computation of the mixed volume of a set of 
bilinear equations using the available software tools pOj can be very demanding. In consequence, the 



exact number of IA solutions is only known for some particular cases [12|, \2l\. 

The main idea that allows us to count efficiently the number of IA solution for single-beam MIMO 
networks is that, as first discussed in (25| , for single-beam MIMO networks the mixed volume does not 
change if we consider rank-one MIMO channels instead of full-rank channels. The proof of this fact is 
straightforward by taking into account that the same monomials are present in both systems of equations 
and, thus, the Newton polytopes that support each equation are identical in both cases. Therefore, for the 
purpose of counting the number of alignment solutions in single-beam feasible systems, we can simplify 
our problem by considering rank-one channels without loss of generality. Assuming rank-one MIMO 
channels, = fkigti, the set of alignment equations ([3]) can be rewritten as 

u|^^3 = o, (M)e*, (20) 

L{u k ) L(vi) 

where now v\ and Uk are column vectors representing the precoders and decoders for the particular case 
of dk = 1 Vfc. We notice that there are K(K — 1) equations, each one being the product of two linear 
factors L(uk) and L(vi) in the entries of u& and v\, respectively, as indicated in ( |20] >. Finding a solution 
to this system reduces to choose from every equation exactly one factor and force it to be zero, i.e., either 
u ]~ fki = or gl t vi = 0. Now the questions is how many different solutions exist for such a system. 

As a first approach, one may think that the total number of solutions would be 2 K ( K ~ 1 ^ since we 
have 2 different choices for each of the K{K — 1) different equations. Actually, this number matches the 
number of solutions given by the classical Bezout's Theorem, which states that the number of solutions 
of a system of polynomials is upper bounded by the product of the degrees of those polynomials. This 
number, however, is just an upper bound (typically very loose) on the number of solutions because the 
equations are coupled through and vi. 

A tighter bound would be obtained by considering that we can design vi (of size Mi x 1)) to lie in, 
at most, Mi — 1 non-intersecting nullspaces. In other words, for a given /, L(vi) = g^vi = can be 
satisfied for a maximum of M; — 1 values of k. This observation would allow us to upper bound the 
number of solutions by or equivalently, (^(7v7-i))' Although this bound is much tighter 

than 2 K<yK ~ l "> for small values of K, they are asymptotically equivalent since the rate of growth with K 
of the latter is also exponential. 

'This is not true though for multibeam cases because in this case the genericity of the system of equations is lost. 
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Due to the fact that all the equations in the system are strongly coupled and they cannot be solved 
independently, the combinatorics of finding the exact number of solutions is much more complicated than 
this last approach and forces us to design a counting routine. In order to explain how this computational 
routine works, we will use the (2 x 3, 1)(3 x 2, 1)(2 x 4, 1)(2 x 2, 1) system as an example. The proposed 
routine proceeds as follows: 

1) We start from a K x K table. Each cell in the table corresponds to a link of the interference 
channel. Cells in the main diagonal represent direct links and they are ruled out since they do not 
play any role in the IA problem. All other cells correspond to interfering links. The table for the 
(2 x 3, 1)(3 x 2, 1)(2 x 4, 1)(2 x 2, 1) system (or any 4-user system) would be as follows. 



Transmitter (I) 





1 


2 


3 


4 


1 










2 










3 






it 




4 











2) We will now fill the cells according to some rules. The value in the cell (k, I) indicates how 
the equation corresponding to the (k,l) link has been satisfied. If it has been satisfied by forcing 
L(vi) = it will contain a one. Otherwise, it will contain a zero, meaning that it has been satisfied 
by setting L{uk) = 0. We recall that given I, L(v{) = can be satisfied for a maximum of M; — 1 
values of k and given k, L(uk) = can be satisfied for a maximum of — 1 values of I. We also 
recall that when s = 0, Y,i( M i ~ l ) + 12k( N k ~ 1) = K{K - 1). Thus, for any valid solution, 
the Z-th column of the table must contain exactly M; — 1 ones whereas the /c-th row must contain 
exactly iVjt — 1 zeros. On the other hand, all cells must contain either a zero or a one. 

As an example, for the (2 x 3, 1)(3 x 2, 1)(2 x 4, 1)(2 x 2, 1) system we just have to find how 
many 4x4 matrices exist with exactly (1, 2, 1, 1) ones in columns 1, . . . , 4 and (2, 1, 3, 1) zeros 
in rows 1, . . . , 4, respectively (not counting those in the main diagonal). It can be seen that there 
are only two possibilities which are shown in Figure [T] 

3) The approach to fill the table for an arbitrary single-beam network is a recursive tree search 
approach, commonly known as backtracking procedure (26} which is widely used to solve combi- 
natorial enumeration problems. We first start with an all-zeros table and try to build up our solution 
cell by cell, filling it with ones, starting from the upper left corner; first right, then bottom. We can 
keep track of the approaches that we explored so far by maintaining a backtracking tree whose root 
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TX antennas (Mi) 
2 3 2 2 



Ones per column 
(M t - 1) 




Solution 1 Solution 2 



Fig. 1. Representation of the two valid solutions for the system (2 x 3, 1)(3 x 2, 1)(2 x 4, 1)(2 x 2, 1). 



is the all-zeros board and where each level corresponds to the number of ones we have placed so 
far. Figure [2] shows the backtracking tree for our example system which was constructed according 
to Algorithm [3] 

In general this procedure is much more efficient than resorting to general software packages to 
compute the mixed volume since it exploits the specific structure of the IA bilinear equations. As 
a matter of fact, our enumeration method for the (3 x 5, l) 7 scenario takes less than 6 minutes to 
count the 357435 solutions, whereas state-of-the-art software for computing mixed volumes could 
take several days to obtain the same solution without additional information regarding the structure 
of the equations. 

1) Connections with other combinatorial and graph theory problems: For the particular case of 
symmetric (M x N, 1) K scenarios, the IA solution counting problem can be restated as several well- 
studied combinatorial and graph theory problems. According to Step 2 above, every column of valid 
solution contains M — 1 ones whereas every row has N — 1 zeros. We recall that s = or, equivalently, 
M + N = K + l, which allows us rephrase Step 2 in a much simpler way. That is, the number of distinct 
IA solutions matches the number of K x K binary matrices with exactly M — 1 ones per row and column 
which are not in the main diagonal. 

When this matrix is seen as the adjacency matrix of a graph (or the biadjacency matrix of a bipartite 
graph) some connections to graph theory problems arise. Most of these problems have been of historical 
interest and hence a lot of research has been done on them. It is natural, then, to find out that the number 
of solutions for some scenarios have already been computed in this field. We mention a few of them in 
the following. 



February 7, 2013 



DRAFT 



17 




Fig. 2. Backtracking tree for the system (2 x 3, 1)(3 x 2, 1)(2 x 4, 1)(2 x 2, 1). The two tables at the lower level correspond 
to the two valid solutions. 



The number of solutions for (2 x (K — i),l) K scenarios is given by the number of derangements 
(permutations of K elements with no fixed points), rencontres numbers or subfactorial. It is also 



the number of labeled 1-regular digraphs with K nodes. Interestingly, as found in |27|, [28 p.195], 
they are equal to 




The number of solutions for (3 x (K — 2), 1) systems matches the number of labeled 2-regular 



digraphs with K nodes. In this case, a closed-form expression is also available |27| 



A (-l) k +j- s Kl(K - k)\(2K - k-2j - s)\ 

Z^Z^Z^ s U k _ s \Wf K _ k _ jy\2j\ 2 2K-2k-j ■ 
k=0 s=0 j=0 y ' vv 

In general, for [M x (K — M + 1), 1) scenarios, closed-form solutions do not exist and most of 
them have not even been studied. It is clear that this problem matches that of counting the number 
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Algorithm 3: Backtracking procedure for counting the number of IA solutions for arbitrary single- 
beam scenarios Yl^ =1 (M^ x Nk,l). 
Input: Number of antennas, {M&} and {N^}; and users, K. 

Output: Number of solutions, 5. 

5 = //No solutions found yet 

table = // Empty table to fill with Is 

row = 0, col = // Row and column indexes 

5 = backtrack {table, row, col, 5) 

function 5 = backtrack (table, row, col, S) 
if table is a valid solution then 

5 = 5 + 1 // Valid solution found 

else 

foreach (row, col) in get_candidates (table,row,col) do 

table(row,col) = 1 // Fill the cell with a 1 

backtrack {table, row, col, 5) // Recursive call 

table(row, col) = // Remove the 1 

return 5 

function ((crowi, ccoh), . . . , (crow^, ccoljy)) = get_candidates (table,row,col) 
return list of candidate cells to store the next 1 



of labeled (M — l)-regular digraphs with K nodes but, as far as we know, no closed-form solution 



has been found yet. Further details can be found in Section IV 



IV. Numerical experiments 

In this section we present some results obtained be means of the integral formulae in Theorem [2] 
(for arbitrary interference channels) and Theorem [3] (for square symmetric interference channels). We 
first evaluate the accuracy provided by the approximation of the integrals by Monte Carlo methods. To 



this end, we focus initially on single-beam systems, for which the procedure described in Section III-D 
allows us to efficiently obtain the exact number of IA solutions for a given scenario. The true number 
of solutions can thus be used as a benchmark to assess the accuracy of the approximation. 
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Tables [IJ and [IT] show the number of solutions given by the exact and the approximate procedures, 
respectively. To simplify the analysis, we have considered (M x (K — M + 1), \) K symmetric single- 
beam networks for increasing values of M and K. As shown in Section III-D1| counting IA solutions 
for this scenario is equivalent to the well-studied graph theory problem of counting labeled (M — Ir- 
regular digraphs with K nodes. Thus, additional terms and further information can be retrieved from 
integer sequences databases such as (27j from its corresponding A-number given in the last row of Table 
[TJ Percentages in Table [IT] represent the upper bound for the relative error, 2e ■ 100, obtained in each 



scenario (see Section III-Cl. 





M = 2 


M = 3 


M = 4 


(2x(jr-i),if 


(3x {K-2),1) K 




K = 2 


i 






K = 3 


2 


1 




K = 4 


9 


9 


1 


K = 5 


44 


216 


44 


K = 6 


265 


7 570 


7 570 


K = 7 


1854 


357 435 


1 975 560 


K = 8 


14 833 


22040 361 


749 649 145 


K > 8 


A000166 


A007107 


A007105 



TABLE I 

Exact number of IA solutions for several symmetric single-beam scenarios, (M x (K — M + 1), 1) 



Figure [3] compares the results in Tables |TJ and [TIJ in a logarithmic scale, to illustrate how the number 
of solutions increase with K for different values of M. The exact number of solutions are also depicted 
to further stress the accuracy of the approximation. 

Now we move to multi-beam scenarios, for which the exact number of solutions is only know for a 
shows the results obtained for some instances of the (M x (2K — M + 2), 2) K 



few scenarios. Table 



III 



network. These results have been obtained using the integral formula in Theorem [2] except the square 
cases (M = N), for which we used the expression in Theorem [3] For instance, we can mention that the 
system (5x5, 2) 4 has, with a high confidence level, about 3700 different solutions. As it can be observed, 
the estimate of the integral formula in Theorem [3] converges much faster than that of Theorem [2j thus 
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M = 2 


M = 3 


M = 4 


(2x (K-1),1) K 


(3x (K-2),1) K 


(4x (K~3),1) K 


K = 


2 


1 


- 


- 


K = 


3 


« 2 (1.0 %) 


1 


- 


K = 


4 


« 9 (1.6 %) 


» 9 (1.6 %) 


1 


K = 


5 


w 44 (2.6 %) 


w 216 (1.5 %) 


w 44 (2.6 %) 


K = 


6 


« 266 (3.3 %) 


w 7 291 (5.5 %) 


w 7 291 (5.5 %) 


K = 


7 


« 1 868 (9.6 %) 


« 361 762 (8.7 %) 


w 1936679 (7.0 %) 


K = 


8 


w 13 144 (20.6 %) 


w 22419610 (11.3 %) 


w 739 668 504 (14.1 %) 



TABLE II 

Approximate number of IA solutions for several symmetric single-beam scenarios, (M x (K — M + 1), 1) 




Fig. 3. Comparison of exact and approximate values for some single-beam scenarios, (M x (K — M + 1), 1) 



allowing us to get smaller relative errors. For the sake of completeness, Table IV shows the approximate 
number of solutions for some additional square symmetric multi-beam scenarios. For some of them the 
exact number of solutions was already known, as indicated in the table. For others (those indicated as 
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N/A in the table) the exact number of solutions was unknown. 



M = 3 M = 4 M = 5 M = 6 



(3x(2^-l),2) K (4 x {2K-2),2) K (5 x (2if - 3), 2) K (6 x (2if - 4), 2) K 
if = 2 1 

if = 3 1 6 1 1 

if = 4 ^9 (5.8 %) w 973 (7.0 %) w 3 700 (0.1 %) « 973 (7.0 %) 

if = 5 ^ 223 (14.8 %) w 530725 (11.3 %) ^ 72581239 (17.8 %) ^ 387 682 648 (0.7 %) 

TABLE III 

Approximate number of IA solutions for several symmetric 2-beam scenarios, (M x (2if — M + 2), 2) K . 



K 


d 


Scenario 


Exact 


Ref. 


Approximate 


3 


1 


(2 x 2, l) 3 


2 




1 




w 2 (0.9 %) 


3 


2 


(4 x 4,2) 3 


6 




1 




w 6 (0.9 %) 


3 


3 


(6 x 6,3) 3 


20 




1 




w 20 (1.4 %) 


4 


2 


(5 x 5,2) 4 


N/A 


N/A 


« 3 700 (0.1 %) 


4 


4 


(10 x 10, 4) 4 


N/A 


N/A 


w 13 887464 893 004 (6.8 %) 


5 


1 


(3 x 3, l) 5 


216 


§ 


w 216 (0.6 %) 


5 


2 


(6 x 6,2) 5 


N/A 


N/A 


w 387 724347 (0.7 %) 



TABLE IV 

Approximate number of IA solutions for selected square symmetric scenarios, (^^d x ^^d, d) 



Although these results have mainly a theoretical interest, they might also have some important practical 
implications. For instance, knowing the rate of increase of the number of solutions with K could have 
interest to analyze the asymptotic performance of linear IA, as discussed in [21]. Also, for moderate- 
size networks for which the total number of solutions is not very high, the results of this paper also 
open the possibility to provide a systematic way to compute all (or practically all) interference alignment 
solutions for a channel realization. Although all IA solutions are asymptotically equivalent, their sum-rate 



performance in low or moderate SNRs behavior may differ significantly [21], [29|. The main idea here 
is that if we are able to obtain all or almost all IA solutions for a particular channel realization, we can 
get all or almost all IA solutions for any other channel realization by using a homotopy-continuation 
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based method such as that described in {25| . This idea is illustrated in Figure |4j which shows in grey the 
sum-rate curves of 973 different solutions for the (4 x 6, 2) 4 network. The maximum sum-rate solution 
is plotted in a thicker solid line, while the average sum-rate of all solutions is represented with a dashed 
line. The relative performance improvement provided by the maximum sum-rate solution over the average 
is always above 10 % for SNR values below 40 dB, and is more than 20 % for SNR=20 dB. We note that 
this improvement is comparable to the one provided by sum-rate optimization algorithms which take into 
account additional information in the optimization procedure such as direct channels and noise variance. 



CO 

+-> 



3 



120 



100 



W 80 



60 



40 



20 - 







All solutions 
Average sum-rate 
Maximum sum-rate solution 




5 10 15 20 25 30 35 40 

SNR [dB] 



Fig. 4. Comparison of the sum rate achieved by approximately different solutions for the system (4 x 6, 2) 4 



V. Conclusion 

In this paper we have provided two integral formulae to compute the finite number of IA solutions 
in tightly feasible problems, including multi-beam (dj. > 1) networks. The first one can be applied to 
arbitrary if -user channels, whereas the second one solves the symmetric square case. Both integrals can 
be estimated by means of Monte Carlo methods. For single-beam networks, it is possible to obtain the 
exact number of solution resorting to more classic results on algebraic geometry. 
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Appendix A 
Mathematical preliminaries 

To facilitate reading, in this section we recall the mathematical results used in this paper. Firstly, we 
provide a short review on mappings between Riemannian manifolds and the main mathematical result 
used to derive the number of IA solutions, which is the Coarea formula. Secondly, we review the volume 
of the complex Stiefel and Grassmanian manifolds and the volume of the unitary group, which are also 
used throughout the paper. 



A. Riemannian manifolds and the Coarea formula 

The following result is immediate from [30) Th. 9.23]. 

Theorem 4: Let X be a compact, embedded, (real) codimension c submanifold of the Riemannian 
manifold Y. Then, for sufficiently small e > 0, 

Vol(y G Y : d(y, X) < e) = Vol{X)Vol(r G R c : ||r|| < l)e c + 0(e c+1 ). 

Here, Vol(X) is the volume of X w.r.t. its natural Riemannian structure inherited from that of Y . 
One of our main tools is the so-called Coarea Formula. The most general version we know may be found 



in 1 31 1, but for our purposes a smooth version as used in |32 p. 241] or |33| suffices. We first need a 
definition. 

Definition A.l: Let X and Y be Riemannian manifolds, and let ip : X — > Y be a C 1 surjective map. 
Let k = dim(y) be the real dimension of Y. For every point x G X such that the differential mapping 
Dip(x) is surjective, let vf, . . . , v% be an orthogonal basis of Ker(Dip(x)) ± . Then, we define the Normal 
Jacobian of <p at x, NJp(x), as the volume in the tangent space T^r x \Y of the parallelepiped spanned 
by Dip(x)(vf), . . . , Dip(x)(vf). In the case that Dp>{x) is not surjective, we define NJ<p(x) = 0. 

Theorem 5 ( Coarea formula ): Let X, Y be two Riemannian manifolds of respective dimensions k\ > 
&2- Let ip : X — > Y be a C°° surjective map, such that the differential mapping Dtp(x) is surjective 
for almost all i£l Let tp : X — ► R be an integrable mapping. Then, the following equality holds: 



4>(x)NJip(x) dX = / if>(x) dx dy. (21) 

xeX J yeY J xGip- 1 (y) 



Note that from the Preimage Theorem and Sard's Theorem (see [34 Ch. 1]), the set ip is a manifold 

is 



of dimension equal to dim(X) — dim(Y) for almost every y G Y. Thus, the inner integral of (21 
well defined as an integral in a manifold. Moreover, if dim(X) = dim(y) then p~ l (y) is a finite set for 
almost every y, and then the inner integral is just a sum with x G p~ 1 (y). 
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The following result, which follows from the Coarea formula, is [32 p. 243, Th. 5]. 

Theorem 6: Let X, Y and V C X x Y be smooth Riemannian manifolds, with dim(V) = dim(X) and 
y compact. Assume that 7r2 : V — > Y is regular (i.e. Dit2 is everywhere surjective) and that Dtt\{x, y) is 
surjective for every (x, y) G V out of some zero measure set. Then, for every open set U C <Y contained 
in some compact set if C 

/ BK 1 W)^= / / DET(x,yy l dxdy, (22) 

JxeC/ ^yey JxeU:(x,y)eV 

where DET(x,y) = det(DG Xty (x)DG Xiy (x)*) and G Xj2/ is the (locally defined) implicit function of tt\ 
near x = 7Ti(x,y). That is, close to (x,y) the sets V and {(x, G Xi y(x))} coincide. 

Corollary 1: In addition to the hypotheses of Theorem [6} assume that there exists yo G Y such that 
for every y G Y there exists an isometry ^ : Y — >• Y with ip y (y) = yo and an associated isometry 
X y : X ->• X such that x y (C7) = U and (x y x <p y )(V) = V. Then, 

/ ${^ 1 {x))dx = Vol{Y) f DET(x, y)^ 1 dx. 

Proof: Let y G Y and let (f y ,Xy as i n the hypotheses. Then, consider the mapping 

Xv l{sel7:(a:,y)eV} : {x G (7 : (x,y) G V} : -)• {x G *7 : (x, y) G V} 

X H- 

which is the restriction of an isometry, hence an isometry. Let G x be the local inverse of iri close to 
x G X. The change of variables formula then implies: 



/ DET(x,y)- 1 dx= f DET( Xy 1 (x),y)- 1 dx. 

J x£U:(x,y)eV J x£U:(x,y )eV 



(23) 



'xGU:(x,y)eV JxeU:(x,y )(: 

Note that the following diagram is commutative: 

TTl it G X) y I 7Tl 

x x 

Thus, the mapping (x^ 1 x ip~ 1 ) o G XjVo o x y is a local inverse of ttl near (x~ 1 (x),y), that is 

Gx-'W,!, = (Xy 1 x V^ 1 ) Ga,, to o 
and the composition rule for the derivative gives: 

^x^*),^ 1 ^)) = X ^y l )^yo^))DG x ,y {x)Dxy Hxy(x))- 

Now, Xy, and Xy x Vy 318 isometries of their respective spaces. Thus, we conclude: 

det(DG x - 1(x) (xy \x))DG x - 1{x)jy ( Xy 1 (*))*) = det(DG x , yo (x)DG x , yo (x)*), 
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that is DET(xZ {x),y) = DET(x,y ). Then, (23 1 reads 



/ DET(x, y)^ 1 dx = / DET(x, yo)' 1 dx. 

Jx£U:(x,y)eV Jx£U:(x,y )eV 



That is, the inner integral in the right-hand side term (22) is constant. The corollary follows 



B. The volume of classical spaces 

Some helpful formulas are collected here: 



27T a 

(cf. pOl p. 248]) Vol(S(C a )) = Vol(S(R 2a )) = —— (24) 

r a 



is the volumen of the complex sphere of dimension a. 



(cf. (35} p. 54]) Vol{U a ) = ) ' (25) 

r(i)-.-r(a) 

is the volumen of the unitary group of dimension a. Note that, as pointed out in (35] p. 55] there are 
other conventions for the volume of unitary groups. Our choice here is the only one possible for Theorem 
|4] to hold: the volume of U a is the one corresponding to its Riemannian metric inherited from the natural 
Frobenius metric in M. a {C). 

We finally recall the volume of the complex Grassmannian. Let 1 < a < b; then, 

Voi(G a , b ) = n^ r(2)-r(6) • (26) 

Appendix B 
Proof of Theorem[T] 

We will apply Corollary [T] to the double fibration given by ( [T0| ). In the notations of Corollary [TJ we 
consider X = T~L, Y = S, V the solution variety and 



^N k -d k J' \0AI k -d k/ 

Given any other element y = (Uk,Vk) E S, let and Qk be unitary matrices of respective sizes iV& 
and Mk such that 

U k = Pk( n Idk ), V k = Q k ( li 



°N k -dJ " \®M k -d, 

Then consider the mapping 

Vy(Uk,V k ) = (P k *U k ,Q* k V k ), 
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which is an isometry of S and satisfies (f y (y) = yo as demanded by Corollary [I] We moreover have the 
associated mapping \y '■ H — > H given by 

Xy{(Hkl)(k,l)<=$) = (Pk H klQl)(k,l)e$ 

which is an isometry of %. Moreover, Xj/C^e) = He and Xy x fyO?) = V- We can thus apply Corollary 
[T] which yields 

[ ^ 1 (x))dx = vol(S) [ (let(DG(H)DG(H)*)- 1 dH, (27) 

where % e is any open subset of % and G is the local inverse of 7Ti close to H at (H,yo). We now 
compute det(DG(H)DG(H)*)~ 1 . From the definition of G we have 

T {H ,u,v)V = {(H, DG(H)H) : H € T H H}. 

On the other hand, from the defining equations ^ and considering H G Hi and H £ Th% as block 
matrices 

. O^xdj A\ < / Ru A kl 

B C J \ Bki Cki 



we can identify 

T iH, yo) v = { (&, ( ° ) = ^ + + Auv t = o 

T/):(tf,F)=-^i? feZ )} 

Hence, 

DG{H)H = -^{R kl ). 
A straight-forward computation shows that: 

DG(H)*(U,V) = (-U^*^ H *(U,V)Vnk^- 
Thus, writing ^ = we have: 

DG{H)DG(H)*(U,V) = * _1 * _ *(J7, F). 
Therefore, (DG(H)DG(H)*)- 1 = and 

detiDGi^DGiH)*)' 1 = det(tf*tf) = det (***). 



From this last equality and (27) we have: 

f %fc\H))dH = Vol{S) [ det(***)diJ. 
Theorem [T] follows dividing both sides of this equation by Vol('H e ) and using the fact that for every 



choice of H out of a zero measure set, the number of elements in ir 1 1 (H) is constant (see [17 Th. 2]). 
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Appendix C 
Proof of Theorem[2] 

Let % e be the product for (k, I) € <5 of the sets 

{H kl : d(H kl ,{R G ^jv fcXMi (C) : ||i2|| F = 1}) < e}. 
From Theorem |4j each of these sets have volume equal to 

2Vol({R e M NkXMl (C) : \\R\\ F = l})e + 0(e 2 ) = r(JVfcMf) e + 0(e 2 ) 



Thus, using ( 24 1, 



(fc,Qe* 

On the other hand, consider the smooth mapping 

/: Hi — > Hi C\X[( kl ^{Hki : \\Hki\\F = ^} 

and apply Theorem [5] to get 

f det(* H *^) dff = / / det(%*l)NJ/(H) dfciH", 

JHeHjOH, J HeH I rt[{^ l y^{HM--\\Hu\\F=l} Jt=(t kl )e[-e,e]»W 

where fljy = H)~i(l + tki)- Note that the function inside the inner integral is smooth and hence for any 
H eHiD ri(jfc,f) 6 ${#W : \\ h m\\f = 1} we have 

det(*£^)NJ/CH) = det(*^)NJ/(H) + 0(e). 

We have thus proved (using « for equalities up to 0(e)): 

/ det(*H*ff)dff« / (2e) tt( * ) det(^^)NJ/(i7)^, 

JHeHtHH, J HeHinU (kl i )e AH M :\\H k i\\ F =l} 

It is very easy to see that NJf(H) = 1 if H = (H^) with ||.ffjy||jr = 1. Thus, we have 

/ det(* ff *^)dH"Rj (2e) s ^ f dei(^ H ^f* H ) dH. 

JHeHrnn, JHeninU (k:l) ^{HM--\\H h i\\ P =i} 

From Theorem [T] and taking limits we then have that for almost every Hq G %, 

t(^\H )) = C-f deb(m*)dH, (28) 

JHeHirin ihm *{H u i\\H u \\r=i} 

where 

2^Vol (5£«inn ( Hw{4 : \\H kl \\ F = l) 
ll(M)e* r(/v fc M,) 
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Now, %i n \\^ k i)£${Hki : \\Hki\\F = 1 is a product of spheres and thus from (24) 

Finally, S = (IljfceT^t.-'V*) x (iLeT ^d^Mi) is a product of complex Grassmannians, and its volume 
is thus the product of the respective volumes, given in (f26|). That is, 



( j lii r(2)-r(iv fe ) J 



<h(.M,-d t ) 



r(2)...r(d { )-r(2)..-r(M,-d l ) 



VtT r(2)..-r(M,) y 

Putting these computations together, we get the value of C claimed in Theorem [2] 

Appendix D 
Proof of TheoremO 

The proof of this theorem is quite long and nontrivial. We will apply Theorem [j] to the sets 

H e = {(H kl ) : d(H kl ,U Nk ) < e,(k,l) G (29) 



Then, because ( 17 1 holds for every e, one can take limits and conclude that for almost every Hq G 



= M«i det(w ., « ( 3o, 

e^O Vol(U e ) JHeUiOfU 

The claim of Theorem [3] will follow from the (difficult) computation of that limit. We organize the proof 
in several subsections. 



A. Unitary matrices with some zeros 

In this section we study the set of unitary matrices of size N > 2d which have a principal dxd minor 
equal to 0, and the set of closeby matrices. For simplicity of the exposition, the notations of this section 
are inspired in, but different from, the notations of the rest of the paper. Let 

T=T N>d =^H=(® d * d A j c M Nx n(C). 

Note that T is a vector space of complex dimension iV 2 — d 2 . Our three main results are: 
Proposition 1: The set Un n T is a manifold of codimension iV 2 inside T ■ Moreover, 

tN-d) 



Vol{U N n T) 



Vol(U N . d f 

Vol(U N -2d) ' 
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Proposition 2: The following equality holds: 

Vol(H eT:d(H,U N )<e) , 



lim 

ei-5-0 



C N 2 



T Vol(U N n T)Vol{x e 



»JV 2 



: \\x\\ < I). 



Proposition 3: Let f : T -> M be a smooth mapping defined in the T and such that ^(H) depends 
only on the A and B part of H, but not on the part C. Denote ^(H) = fy(A,B). Then, 

I H eT:d(H M <e^( H ) dH 



!™ Vol{H £ M N (C) : d(H,U N ) < e) 



2 d2 Vol{U 



N-d 



)f(A 



Vol (U N ) Vol {U N -2d)J(A' ,B)eu { _ 
1 ) Proof of Proposition [TJ Let 



V(A,B)d(A,B) 



e: U 2 N _ d U N nT 



(U,V) ^ 



Id (> 

u, 



J 



Id o 
V* 



(31) 



where 



J 

We claim that £ is surjective. Indeed, let 

# = 



J d 
Id 
I N . 2d ) 



A 
B C, 



eU N nT. 



From HH* = In we have that A satisfies AA* = Id, i.e. the rows of A can be completed to form a 
unitary basis of C N ~ d . Namely, there exists V G U^-d such that A = (Id 0)V. Similarly, there exists 
U £ U N _ d such that B = U(^). Then, 



H 



h 

u 



1 o i d o N 



Id Rl R'2 

\0 Rs Ra) 

where R satisfies U RV = C. Now, this implies that the matrix 

^ o i d o N 

Id Ri R2 
\0 R 3 R 4 J 



Id 
F 
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is unitary, which forces Ri = 0, R 2 = 0, R3 = and i? 4 unitary. That is 

/n t, n \ 



H 




I d 
Id 

v r 4 J 





( h 0^ 

I d 

y i? 4 y 




that is 



H = £ U,V* 



Id 
#4 

and the surjectivity of £ is proved. Moreover, this describes Un n T as the orbit of J under the action 
in T given by 

'id 




V* 



Then, £Yjv H T is a smooth manifold diffeomorphic to the quotient space 

Uff_ d /Sj, 



where Sj is the stabilizer of J. Now, (U, V) £ Sj if and only if 



fid \ 

C/i c/ 2 

yo J7 3 ^ 4 y 



/ d 
h 

I N .2d 



h 



\ 



Kj* 



V? Va 



4/ 



i d 
id 
001 






N—2d 



which implies Ui = h, U 2 = 0, U 3 = 0, V\ = I d , V 2 = 0, V 3 = and f/ 4 = F 4 . Thus, 



Sj 




: U4 £ UN-2d 



(32) 



Then, 



dim(Wjv n T) = d\m{U 2 N _JSj) = 2 dim(U N _ d ) 2 - dim(Sj) = 2(N - d) 2 - (N - 2d) 2 = N 2 — 2d 2 . 

On the other hand, dim(T) = 2N 2 - 2d 2 and thus 

codim r (^iv n T) = 2N 2 - 2d 2 - (N 2 - 2d 2 ) = N 2 , 

as claimed. We now apply the Coarea formula to £ to compute the volume of Un H T ■ Note that by 
unitary invariance the Normal Jacobian of £ is constant, and so is Vol{^~ l {H)). We can easily compute 

V2 {N - 2d)2 vom 



Vol{C l (H)) = VoliCV)) = Vol(Sj) = 

V H (K32l 



N-2d 
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D£(I N - d ,I N - d ){U,V) 



\ 



For the Normal Jacobian of £, note that 

( V{ -v 2 
iii o u 2 
\-u* 2 v 2 * Ui + vz) 

Thus, D^ijjsi-d, In-cl) preserves the orthogonality of the natural basis of TulA^-d x TyU^-d but for 
the elements such that U4 ^ or V4 ^ 0. We then conclude that iVJ(^)(Jjv_^, In-<i) = NJ{rj) where 

77: {M G M N - 2d ■ M + M* = 0} 2 -> {M £ M N - 2d - M + M* =0} 

It is a routine task to see that rj*(L) = (L, L*) which implies r}rj*(L) = 2L, that is 

detOqr]*) = 2 dim ^ MeMN - 2d - M+M "= ^ = 2( Ar " 2d ) 2 . 

Hence, NJ(rj) = \J det (rjr/*) = \f2^ N 2< ^ . As we have pointed out above, the value of the Normal 
Jacobian of £ is constant. Thus, for every U, V, 

NJ(£)(U,V) = NJ( V ) = V2 (N ' 2dr . 



The Coarea formula applied to £ then yields: 



Vol{U 



N-d) 



(v,v)mi 



1 d(U, V) 



Heu N nT 



NJ(0 



dH = Vol(U N n T)Vol{U N -2d)- 



The value of VoI(Un H T) is thus as claimed in Proposition [I] 
2) Some notations: Given a matrix of the form 

/n . n\ 



H 



a 

a Ci C 2 

\o c 3 c 4 y 



(33) 



(a and a are d x d diagonal matrices with real positive ordered entries) we denote by H the associated 
matrix 



H 



\ 



a Ci C 2 
cr 

where Uq is some unitary matrix which minimizes the distance from C4 to U^-id- Note that 

\ 



H 



^0/0 
7 
V° ^/ 



H, 
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and hence 

d(H,U N ) = d{H,U N ). 



We also let 



rr(n\ IU ri|2, \\„ ri|2, nr. tt ii2 , ll^ill 2 + IIC2II 2 + \\C3W" 
±i(H) = \\a - ld\\ +||cr — id || +HC4 — l/o|| H 



Note that 



2 

r 2 (ff) = \\ a - i d \\ 2 + Ik - / d |l 2 + l|C7 4 - u f + HdH 2 + IW + IW 



m/r« m ,rrt II # ~ -^AT 1 1 2 d(H,U N ) 2 d(H,U N ) 2 

T 2 (H) > Ti(H) > * > V 2 ^ = V 2 W; (34) 



3) Approximate distance to Un and Un H T: In this section we prove that for small values, 

d(H,U N ) ~ Ti(H) 1 / 2 , d(H,UN nf)« T 2 (F) 1/2 . 

More precisely: 

Proposition 4: For sufficiently small e > 0, if d(H,UN) < e then, 

IdiHM-T^H) 1 ^ <0(e 2 ), 

d(H,U N nT) -^{H) 1 ^ < 0(e 2 ). 

Here, we are writing 0(e 2 ) for some function of the form c(d)e 2 . 
Before proving Proposition [4] we state the following intermediate result. 
Lemma 1: There is an eo > such that \\H — < e < eo implies: 

Ti(i?) 1/2 - 9e 2 < d(H,U N ) < Ti(#) 1/2 + 9e 2 , 

T 2 (#) 1/2 - 30e 2 < d(H,U N r\T) < T 2 (H) 1/2 + 30e 2 . 



Proof: We will use the concept of normal coordinates (see for example [30 p. 14]). Consider the 
exponential mapping in Un, which is given by the matrix exponential 

TfU N = {R G M N (C) : R + R* = 0} -»■ Wjv 

fl _> e« = / + fl + E fea ^, 

which is an isometry from a neighborhood of G TjWat to a neighborhood of I G Wat and defines 
the normal coordinates. Thus, for sufficiently small e\ > there exists eo > such that if £7 G IAn, 
\\U — J|| < eo then there exists a skew-symmetric matrix i? such that 

U = e R , \\R\\=d UN (U,I), ||fl||<ei. 
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Let R G .Mjv(C) be a skew-Hermitian matrix such that 

\\H-e R \\ = d(H,U N ) = S<e, \\R\\ = d UN (e R ,I), \\R\\ < t x . 
Let S = Y, k >2 Rk / kl Then ' e R = I + R + S and 

\\ S \\<^^t<2\\R\\ 2 . 



k>2 



k\ 



If we denote a = \\e R — I\\ = ||-R + S\\ and b = du N (e R ,I) = \\R\\, we have proved that 

b-2b 2 <a<b + 2b 2 . 
Using that b < 1/2 and some arithmetic, this implies 

a + 6a 2 >b, that is \\R\\ < \\e R - I\\ + 6\\e R - I\\ 2 . 

Now, 

\\e R - I\\ < \\e R - H\\ + \\H - I\\ < 2e, 

which implies 

p|| < 2e + 24e 2 < 3e. 

In particular, ||5|| < 9e 2 . We conclude: 

d(H,U N ) = d(H,U N ) = \\H - e R \\ > \\H - (I + R)\\ - \\S\\ > \\H - (I + R)\\ ■ 
We now solve the following elementary minimization problem: 

min \\H-(I + R)\\. 

R:R+R'=0 



Let 



R = 



( R\ i?2 

-R% i?5 Ra 



\-H 3 -Kg 



, R\ + R* — 0, i?5 + i?5 — 0, Rg + i?g — 0. 



R9J 



Then, \\H — (I + R)\\ is minimized when R± = 0, R5 = 0, Rg = and 

R 2 =argmin(\\Ci - R 2 \\ 2 + \\R2f) 
R 3 =argmin{\\C 2 - Rzf + \\Rzf) 
R 6 =argmin(\\U^C 3 - R* 6 \\ 2 + ||i? 6 || 2 ). 
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It is easily seen that the solutions to these problems are: 

Cl . t-, ii 2 , ii t-, ii 2 IIClll 2 



#2 
#3 
i?6 



2 

" 2 ^ 
C* 3 U 



I Ci ~~ -R2 1 1 2 + II R2 1 1 2 



2 ' 



1 " 1 1 2 1 1 2 II ^2 

|C2 — i?3 + i?3 = — — 



\mc 3 -m\\ z + \\R 6 



M,2 . iini 2 licfcll 1 



We have then proved 



min \\H - (I + R)\\ = TxiH) 1 ' 2 , 



and the minimum is reached at 



R 



Hence, 



( Ci/2 C 2 /2 N 

-Cf/2 Wo/ 2 
r C*/2 -C/oq/2 / 



d(H,U N ) > Ti(H) 1 / 2 — 9e 2 , 



(35) 



and the first lower bound claimed in the lemma follows. For the upper bound let R be defined by (35 1 
and note that (following a similar reasoning to the one above) 



d(H,U N ) = d(H,U N ) < \\H-e R \\ < \\H-(I+R)\\+^2 



\R\\ 



k\ 



T 1 (H) 1 / 2 +^2 



k>2 k>2 

12 1 \\n 112 1 \\n 112 s A 



2 

fc! 



k/2 



Now, \\H — Jiv|| < e in particular implies ||Ci || + \\C2W + IIC3II — e an d then we have 

k/2 



d(H,U N ) < T^H) 1 ' 2 + £ V < T^tf) 1 / 2 + 2e 



fc>2 



as wanted. Now, for the second claim of the lemma, the same argument is used but now R is such that 
e R minimizes \\H — e^ll and 



Now, from the equality 



* * 



I + R = e R-S 



and arguing as above we have that 

\\R2W < \\S\\ < 9e 2 , which implies \\R - R\\ = y/2\\R 2 \\ < 20e 2 , 
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where we denote by R the matrix resulting from letting R2 = 0. Thus, 

lll-fiT-e^H - \\H-(I + R)\\ \ <\\\H-I-R\\- \\H - I - R\\ | + ||5|| < \\R - R\\ + 9e 2 < 30e 2 . 
We have then proved 



d (H,Un H 7~) — min \\H - (I + R)\\ 
R:R+R'=0,R 2 =0 



< 30e 2 , 



and as before we can easily see that the minimum is reached when R\ = 0, R2 = 0, R$ = 0, Rg = 0, 
R 3 = C 2 /2 and R 6 = C%U /2 which proves that 

min \\H-(I + R)\\=T 2 (H). 
R:R+R*=0,R 2 =0 

The lemma is now proved. ■ 
Proof of Proposition [?] 

Let E be a matrix such that ||i?|| < e and H = U + E for some unitary matrix U. Then, 

\\HH* - III = \\UU* + UE* + EU* + EE* - III = II [IE* + EU* + SS* II < 2e + e 2 < 4e. 



On the other hand, 



(a 2 -I oC\ aCl \ 



HH* - I 



da X X 
\ C 3 a X C 3 CI + C 4 Ct-IJ 
where the entries X are terms which we do not need to compute. In particular, we have ||Ci<t|| < 4e and 

||cj 2 -I|| < 4e, (36) 

which implies ||o- _2 || = 1 1 cr 2 — I + I|| < V~d + 4e and hence 

||d|| = Wdaa^W < HCictIIIIct- 1 !! < Aey/y/d + Ae < iVde. 

A similar argument works for C3 as well, and using a symmetric argument for H*H we get the same 



bound for C2 and an equivalent bound for a to that of (36 1. Summarizing these bounds, we have: 

||Ci|| 2 + ||C 2 || 2 + ||C 3 || 2 <48de 2 (37) 

Moreover, we also have 

\\C4Cl - I|| < HC3C3 || + HC3C3* + C4C4* - I|| < 16de 2 + 4e < 20de, 

which implies 

N—d 

(Pj - l) 2 = \\C A Cl - I|| 2 < 400d 2 e 2 

j=0 
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where the /3j are the singular values of C4. In particular, 

N-d N-d N-d 

\\UZC 4 -I N - d \\ 2 = d(C A ,U N - d ) 2 = ^(/3,-l) 2 < ^(/3,-l) 2 (& + l) 2 = ]T(/3 2 -l) 2 < 400d 2 e 2 

j=i j=i j=i 

and we conclude that 



\UoC4-I N - d \\ < 20de. 



(38) 



Using (36 1, (37 1 and (38 1 above we get: 



\\H - I N \\ 2 = \\a - I d \\ 2 + \\a - I d \\ z + \\Ctf + \\C 2 \\ Z + \\C 3 \\ Z + ||C/ *C 4 - I N - d \\ 2 < C(d) z e z , 

where C(d) depends only on d. Let e be small enough for C(d)e to satisfy the hypotheses of Lemma [T] 
The Proposition [4] follows from applying that lemma. 

4) How the sets of closeby matrices to Un and Un H T compare: Our main result in this section is 
the following. 

Proposition 5: Let a > 1. For sufficiently small e > 0, we have: 



2 d Vol \H G T : d(H,U N n T) < -J < 
yoZ(fT eT:d(H,U N ) < e) < 
2 d2 Vol{H G T : d(fl",Wjv n T) < ae) 



Before the proof we state two technical lemmas. 



Lemma 2: Let a, a be as in ([33]). Then, 
/ 



Proof: Let 



C : Ti 



a 

a Ci C 2 

5i(C7) 



< e 



C : T 2 



<r 

a Ci C 2 

\o c 3 c 4 y 



< e 



/ 



A 
B C, 



1,2, 



where i = (d 0) and 1? T = (a 0). The claim of the lemma is that 

Vol(C : S'i(C) < e) = 2 d Vo/(C : 5 2 (C) < e). 
Indeed, consider the mapping 



C\ C 2 
C3 c 4 



V2d C 2 
C3 C 4 
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which has Jacobian equal to \ff' d = 2 d2 . The change of variables theorem yields: 

2 d2 Vol{C : 5i((^(C)) < e) = Vol{tp(C) : S^C)) £ e) = Vol(C : Si(C) £ e). 

The lemma follows from the fact that Si(<p(C)) = S 2 (C). ■ 
Lemma 3: Let a > 1 and let A, B be complex matrices of respective sizes dx (N — d) and (N — d) x d. 
Then, for sufficiently small e > we have 

2 d2 Vol | C : d | | ° A ^ Mn n T^J £ ^ 

FoZ | C : d | | ° | M N | ■ | 

B C, 

2^ Vol | C : d | I ° | .//v nT 1 £ n< 
\B C, 

Proof: Let Ua,Va,Ub, Vb be such that 

= Ua(o~ 0)V£, B = U B (^jV B 
are singular value decompositions of A and B respectively. Then, 




A\ IV B 
B CI \ Va, 





= FoZ ( 







I V,/ If: J I | U I ,Wjv I £ e I = FoZ [ C : d \ | ° I .//\ I £ < 



where the last inequality follows from unitary invariance of the volume. Let H be as in ([33]). From 
Proposition |4| we conclude: 

Vol (C : d{H,UN) < e) < ^oZ(C : T^H) 1 ' 2 < e + c((i)e 2 ) = 

Vol(C : Ti(H) £ (e + c(d)e 2 ) 2 ) = rn 2 d Vo/(C : T 2 (#) £ (e + c(d)e 2 ) 2 ). 

Lemma|2] 



From (34 1, for sufficiently small e > 0, T 2 (#) £ (e + c(cf)e 2 ) 2 implies d(H,U N ) is as small as wanted. 
Hence, from Proposition |4j for sufficiently small e > we have 

Vol(C : T 2 {H) £ (e + c(d)e 2 ) 2 ) = FoZ(C : T 2 {H) 1 ' 2 < e + c(cZ)e 2 ) £ 

Vol (C : d (H, Un n 7") £ e + 2c(d)e 2 ) . 
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In particular, for every a > 1 and for sufficiently small e > we have proved that 

Vol(C : d(H,U N ) < e) < 2 d VoZ (C : d{H,U N nT) < at). 

This proves the upper bound of the lemma. The lower bound is proved with a symmetric argument, using 
the opposite inequalities of Proposition [4] 

■ 

Proof of Proposition [5] Let a > 1. From Fubini's Theorem, 

Vol{H € T : d(H,U N ) < e) = / y /(C : d(H,U N ) < e) £). 

JAEMdx(_N-d) (C),BeX ( jv-d)xd(C) 

From Lemma [3j for sufficiently small e > this is at most 

[ 2 d2 Vol(C :d(H,U N nT)<ae)d(A,B). 

(C) 

Again from Fubini's Theorem, this last equals 

2 d2 Vol{H : d(H,U N n T) < ae), 

proving the upper bound of the proposition. The lower bound follows from a symmetrical argument. 

5) Proof of Proposition^ Let a > 1. From Proposition [5] we have 

, Vol(HeT:d(H,U N )<e) # y Vol(H £ T : d(H,U N CiT) < ae) 
lim ^ < 2 lim 



e^o e N2 - e£o e N2 

Note that N 2 is the (real) codimension of Un H T inside T. Thus, from Theorem |4j 

lim ^ = VoliptN n T)a Vol[x G K : < 1). 

We have thus proved that for every a > 1 we have 

Hm Vol(HGT:d(H,U N )<e) ^ p e R *» . || x || < 1} . 

This implies: 

^ Vol{H^T:d{HMN)<e) < p £ ^ ; < 1} 

The reverse inequality is proved the same way using the other inequality of Proposition [5] 
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6) Integrals of functions of the subset of matrices in T which are close to U^: We are now close to 
the proof of Proposition [3j but we still need some preparation. We state two lemmas. 
Lemma 4: Let ^ : T — > [0, oo) be a smooth mapping. Then, 



^ / $>(H) dH = 2 d2 Vol(U N n T)Vol{x G R N ' : ||x|| <l)-f V(U) dU 

JH£T:d(H,u N )<e JUeU K nT 



Proof: For sufficiently small e > 0, given H G T such that d(H,UN) < e, there is a unique 
U G U n T such that the distance d(H,Un T) is minimized (see for example |30| p. 32]). Let ir(H) be 
such U. Moreover, it is a smooth mapping. We thus have 



L 



^dH= [ [ NJir(H)^(H)dHdU. 

JUGU^nT JHeT:d(H.U^)<€.TT(H)=U 



'H&T:d(H,U N )<e JUeU N nT J H£T:d{H,U N )<e,Tr(H)=U 

Now, * is smooth and hence ^>{H) =<S{U) + ^>{H) where \&(H) | < O(e). We thus have 



$>dH= [ V(U) [ NJir(H)dHdU + 0(e) 



'HeT:d(H,U N )<e JU&U N nT J H&T:d(HJJ N )<<i,-K{H)=U 

The integral inside this last expression is unitary invariant and thus its value is a constant c e . Moreover, 
the same argument applied to Vl/ = 1 yields 



Vol(H G T : d(H,U N ) < e) = / c e dU. 



That is, 

_ Vol(H eT:d{H,U N ) < e) 
Ce ~ Vol{U N n T) ' 

We have then proved 



y dH= Vol(HeT:d(H,m<<) ! nu)du 
HeT-.d(H,u N )<e Vol(U N D T) Jueu N nr 

Vol(H GT : d(H,U N ) < e)i ^{U)dU. 

JupUmDT 



The lemma follows from Proposition [2j 

Lemma 5: Let $ be a smooth mapping. Then, 

I H er A H W <MH)dH = ^Vol(U N nT) f 
e^oVol(H GM N (C):d(H,U N )<e) Vol{U N ) J u& A N nT 

Proof: From Theorem [4] and using that the codimension of Un in M.^^) is N 2 we know that 

Vol(H G M N {C) : d{H,UN) < e) = Vol{U N )e N2 Vol{x G R N " : \\x\\ < 1)(1 + 0(e)), 

where lim el _>o 0(e) = 0. The lemma now follows from Lemma |4] 
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7) proof of Proposition [£]■ This result is almost inmediate from Lemma [5] and Proposition [T] Let £ 
be the mapping defined in pT) . We have computed the Nromal Jacobian of £ and the volume of the 
preimage of £ in Section D-Al From Theorem [5] 



/ vfr(e(Z7, VO) V) = I ^(H) V0l{ l {H)) dH = Vol(U N - 2d ) [ 9(H) dH 

J(U,V)&u%_ d JHeU N r\T iVJ ? JHeU N nT 

Hence, as *J7 does not depend on C, and writing 9(H) = 9(A,B) (note the abuse of notation), 

9 W dH = i , I 9((I d 0)V*,u( I A) d(U,V). 

H&u N nT Vol(U N -2d) J(u,V)eU 2 N _ d V V u // 

Normalizing we get 

/ 9(H)dH = l 9((I d 0)V*,u( 1 *)) d(U,V). 

JH&u N c\T J{uy)&A 2 N _ d \ \ u // 

Now, generating at random unitary matrices U,V and then taking (I d 0)F*,?7(Q d ) is the same as 

generating at random two elements in Stiefel^_ d ^ xd . The proposition is proved. 

B. Proof of Theorem [J] 



Recall that we have defined % e in (29), and we want to compute the limit (30 1 



,. Vol(HinH e )Vol(S) f , , _ » ,. Vol(S) f , 

lim v - ttztts —4 det(99*) dH = lim — — )-/- / det (99*) dH. 

JHeHrnn, eh ^° Vol(H e ) JH&i,rOi, 



Vol(H e ) JHeHjnHc Vol(H e ) Jhzh^h, 

Now, we use Fubini's theorem to convert the last integral into an iterated integral 



det(99*)dH {krjlr y dH {kuh) , 



where (k±,li), . . . , (k r ,l r ) are the elements of the connectivity set 9, ordered with respect to some 
(irrelevant) criterion. From Proposition |3j the last inner integral satisfies: 



/ det(99*) dH {kr>lr) = 0(e*) + Vol(H G M N (C) : d(H,U N ) < e) 

det(W) d(A,B), 



„. - ) 

' H( kl ^i 1 )£_T,d(H( kr ^ r .),U N )<E 

2 d2 Vol(U N . d ) 2 



Vol(U N )Vol(U N -2d)J(A*,B)GU iN - d)xd 

where 9 is computed for 

'Odxd A 



H {kr,3r) 



Here, 0(e*) is an expression such that 



B Q(N-d)x(N-d) 



lim m = 0. 

e^o Vol(H e M N (C) : d(H,U N ) < e) 
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By repeating the procedure and using Fubini's theorem again to convert the iterated integral into a unique 
multiple integral, we conclude: 

/ det(W)d# = 0(e*) + Vol(H € M N (C) : d(H,U N ) < e) m x 



2 d 'Voi(u N - d y 



7 h 



Vol{U N )Vol{U N -2d) J J(A' M ,B M )eU (N - d)xd ,(k,i)e^ 
where \t is computed for 



if, 



ki 



B k i 0{N-d)x{N-d) 



Here, 0(e*) is an expression such that 



hm m = o 

e^o Vol(H e M N (C) : d(H,U N ) < e)»(*) 
On the other hand, also from Fubini's theorem we have 

Vol(H £ ) = Vol(H G M N (C) : d(H,U N ) < e)^l 

The claim of the Theorem [3] follows. 
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